Search for: All records

Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

 

We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

 

Understanding the temporal variability of plate tectonics is key to unraveling how mantle convection transports heat, and one critical factor for the formation and evolution of plate boundaries is rheological “memory,” that is, the persistence of weak zones. Here, we analyze the impact of such damage memory in global, oceanic‐lithosphere‐only models of visco‐plastic mantle convection. Self‐consistently‐formed weak zones are found to be reactivated in distinct ways, and convection preferentially selects such damaged zones for new plate boundaries. Reactivation of damage zones increases the frequency of plate reorganizations, and hence reduces the dominant periods of surface heat loss. The inheritance of distributed lithospheric damage thus dominates global surface dynamics over any local boundary stabilizing effects of weakening. In nature, progressive generation of weak zones may thus counteract and perhaps overcome any effects of reduced convective vigor throughout planetary cooling, with implications for the frequency of orogeny and convective transport throughout Wilson cycles.

 

Azimuthal variations in receiver function conversions can image lithospheric structural contrasts and anisotropic fabrics that together compose tectonic grain. We apply this method to data from EarthScope Transportable Array in Alaska and additional stations across the northern Cordillera. The best-resolved quantities are the strike and depth of dipping fabric contrasts or interfaces. We find a strong geographic gradient in such anomalies, with large amplitudes extending inboard from the present-day subduction margin, the Aleutian arc, and an influence of flat-slab subduction of the Yakutat microplate north of the Denali fault. An east–west band across interior Alaska shows low- amplitude crustal anomalies. Anomaly amplitudes correlate with structural intensity (density of aligned geological elements), but are the highest in areas of strong Cenozoic deformation, raising the question of an influence of current stress state. Imaged subsurface strikes show alignment with surface structures. We see concentric strikes around arc volcanoes implying dipping magmatic structures and fabric into the middle crust. Regions with present-day weaker deformation show lower anomaly amplitudes but structurally aligned strikes, suggesting pre-Cenozoic fabrics may have been overprinted or otherwise modified. We observe general coherence of the signal across the brittle-plastic transition. Imaged crustal fabrics are aligned with major faults and shear zones, whereas intrafault blocks show imaged strikes both parallel to and at high angles to major block-bounding faults. High-angle strikes are subparallel to neotectonic deformation, seismicity, fault lineaments, and prominent metallogenic belts, possibly due to overprinting and/or co-evolution with fault-parallel fabrics. We suggest that the underlying tectonic grain in the northern Cordillera is broadly distributed rather than strongly localized. Receiver functions thus reveal key information about the nature and continuity of tectonic fabrics at depth and can provide unique insights into the deformation history and distribution of regional strain in complex orogenic belts. 

The spatial distribution of the geochemical domains hosting recycled crust and primordial (high‐³He/⁴He) reservoirs, and how they are linked to mantle convection, are poorly understood. Two continent‐sized seismic anomalies located near the core‐mantle boundary—called the Large Low Shear Wave Velocity Provinces (LLSVPs)—are potential geochemical reservoir hosts. It has been suggested that high‐³He/⁴He hotspots are spatially confined to the LLSVPs, hotspots sampling recycled continental crust are associated with only one of the LLSVPs, and recycled continental crust shows no relationship with latitude. We reevaluate the links between LLSVPs and isotopic signatures of hotspot lavas using improved mantle flow models including plume conduit advection. While most hotspots with the highest‐³He/⁴He can indeed be traced to the LLSVP interiors, at least one high‐³He/⁴He hotspot, Yellowstone, is located outside of the LLSVPs. This suggests high‐³He/⁴He is not geographically confined to the LLSVPs. Instead, a positive correlation between hotspot buoyancy flux and maximum hotspot³He/⁴He suggests that it is plume dynamics (i.e., buoyancy), not geography, which determines whether a dense, deep, and possibly widespread high‐³He/⁴He reservoir is entrained. We also show that plume‐fed EM hotspots (enriched mantle, with low‐¹⁴³Nd/¹⁴⁴Nd), signaling recycled continental crust, are spatially linked to both LLSVPs, and located primarily in the southern hemisphere. Lastly, we confirm that hotspots sampling HIMU (“high‐μ,” or high²³⁸U/²⁰⁴Pb) domains are not spatially limited to the LLSVPs. These findings clarify and advance our understanding of deep mantle reservoir distributions, and we discuss how continental and oceanic crust subduction is consistent with the spatial decoupling of EM and HIMU.

 

Subducted sediments are thought to lubricate the subduction interface and promote faster plate speeds. However, global observations are not clear‐cut on the relationship between the amount of sediments and plate motion. Sediments are also thought to influence slab dip, but variations in subduction geometry depend on multiple factors. Here we use 2D thermomechanical models to explore how sediments can influence subduction dynamics and geometry. We find that thick sediments can lead to slower subduction due to an increase of the megathrust shear stress as the accretionary wedge gets wider, and a decrease in slab pull as buoyant sediments are subducted. Our results also show that larger slab buoyancy and megathrust stress due to thick sediments increase the slab bending radius. This offers a new perspective on the role of sediments, suggesting that sediment buoyancy and wedge geometry also play an important role on large‐scale subduction dynamics.